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Initially produced quark-gluon matter at RHIC and LHC does not have a temperature. A 
quark-gluon plasma has a high temperature. From this quark-gluon matter to the quark- 
gluon plasma is the early thermalization or the rapid creation of temperature. Elastic 
three-parton scattering plays a key role in the process. The temperature originates from 
the two-parton scattering, the three-parton scattering, the four-parton scattering and so 
forth in quark-gluon matter. 
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Temperature is common in daily life, and people believe that temperature exists 
no matter how cold matter is. In gas and liquid temperature reflects the existence 
of thermal states where all particles randomly move and obey isotropic momentum 
distributions. The average momentum that a massless particle possesses is propor¬ 
tional to temperature. When temperature equals zero, the average momentum is 
zero, and vice versa. However, there is one system where the average momentum is 
large but particles do not obey isotropic momentum distributions, i.e., the system 
does not have temperature. The system is quark-gluon matter that was produced in 
initial gold-gold nuclear collisions performed at the Brookhaven National Labora¬ 
tory Relativistic Heavy Ion Collider (RHIC) in the last fourteen years. In the parton 
system the parton distribution functions in the direction of the incoming nucleus 
beam are much larger than those in the direction perpendicular to the beam direc¬ 
tion. The distributions are anisotropic in momentum space, and do not correspond 
to a thermal state. Now we must wonder if this quark-gluon matter always has no 
temperature. The answer is no, and this matter is exactly a place for exploring the 
origin of temperature. 

The reaction gg —?► gg was considered by Shuryak— to cause gluon matter to 
thermalize. He assumed that every gluon suffers from scattering once in the same 
period of time. In order to guarantee validity of the assumption, the cross section 
for the gluon-gluon scattering must be very large (empirically be larger than 40 
mb). This is not achievable! Even though a thermalization time of the order of 0.3 
fm/c was obtained, the assumption makes the thermalization time unreliable. To 
explore the kinetic equilibration of gluon momenta, inelastic scattering gg —^ ggg 
was taken into account in Refs .— gg —>■ gg and gg —^ ggg lead to a gluon-matter 
thermalization time larger than I fm/c. Overpopulation is found to be a key feature 
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of the Glasma (far-from-equilibrium gluon matter), and can strengthen the kinetic 
equilibration. But solutions of a transport equation with the elastic gg scattering 
show that the thermalization time of the Glasma is also larger than 1 fm/o^. 

However, the kinetic equilibration, i.e., the establishment of a thermal state 
takes a time less than 1 fm/c, as concluded from the elliptic flow data of hadrons 
produced in Au+Au collisions at y/s_NN = 130 GeV at RHIG^ and the correspond¬ 
ing explanation of hydrodynamic calculations which assume early thermalization of 
quark-gluon matter and ideal relativistic fluid flow of the quark-gluon plasma^ii. It 
was pointed out by Kolb et al. that quark-gluon matter begins to have a temper¬ 
ature at 0.6 fm/c from the moment when quark-gluon matter is created in initial 
Au-|-Au collisions. The early thermalization thus became the first RHIC puzzle^. 
Even for Pb-|-Pb collisions at ^/snn = 2.76 TeV at the Large Hadron Gollider 
(LHG), a good explanation of the elliptic flow data of hadrons^ in hydrodynamic 
calculations requires that quark-gluon matter begins to have a temperature at 0.6 
fm/o^. Therefore, the early thermalization, i.e., the rapid equilibration also exists 
at LHG. 

To explain the early thermalization or to find the origin of temperature of quark- 
gluon plasma, we think about the prominent feature of initially produced quark- 
gluon matter: the parton number density is as high as 32 fm~^. At such a high num¬ 
ber density the occurrence probability of parton-parton-parton scattering is larger 
than 0 . 2 ^. The feature leads to that the three-parton scattering is important in es¬ 
tablishing thermal states of quark-gluon matter, i.e., creating an initial temperature 
of quark-gluon plasm a^^d^ . For example, gluon-gluon-gluon scattering was discov¬ 
ered to dominate the evolution of gluon matter towards thermal equilibrium—. We 
draw the conclusions from transport equations which include the elastic parton- 
parton scattering— and the elastic scattering of gluon-gluon-gluon—, gluon- 
gluon-quark, gluon-gluon-antiquark—, gluon-quark-quark, gluon-quark-antiquark, 
gluon-antiquark- antiquark—, quark-quark-quark—, quark-quark-antiquark, quark- 
antiquark-antiquark—, and antiquark-antiquark-antiquark. The transport equation 
for gluon matter is 

^ _1 f (Pp2 d^P3 d^P4 .4 

dt ^ 2Ei J (27r)32A2 {2tt)^2E3 (27r)32A4 ^ ’ 

xS'^ipi +P2-P3-P4 )[y I -^SS^gg [/gl/g2(l + /g3)(l + /g4) 
-/g3/g4(l + /gi)(l + /g2)] + gqd Algu^gu P + I A^gd^gd P 

+ I Afgu-J-gu I + I AIgd_>gd I )[/gl/q2(l + /g3)(l ~ /qd) 

-/g3/q4(l + /gl)(l - /q2)]} 

_^ f <Pp2 d?P3 (Ppi d?Pb (PpG 

2Ei J (27r)32A2 (27r)32E3 (27r)32A4 {2'kY 2E^ {2Tr)32Ee 

X (27r)4(5'‘(pi +P2+P3-P4-P5-P6)1%\ Alggg^ggg P [/gl/g2/g3 
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x(l + /g4)(l + /g5)(l + /ge) - /g4/g5/g6(l + /gl)(l + /g2)(l + /gs)] 

H 2 (I -^ggu-s-ggii I I -^ggd^ggd I + I -^ggu-^ggu I + I -^ggd-^ggd I ) 

X [/gl/g2/q3(l + /g4)(l + /g5)(l ~ /qd) ~ /g4/g5/q6(l + /gl)(l + /g2) 

X(1 — /qs)] +5 q[^ I -A^Iguu-^guu P “*" 2 ^^ -^gud-s-gud P + | A^gdu-)-gdu P) 

”^4 ^ -^gdd-J-gdd I + I -^guu-^guu | + | -^gud-fgud I + I -^gdu-s-gdu | 

+ I -^gdd^-gdd P I -^guu-^guu P + 2(1 -^gud-i-gud P + I -^gdu-s-gdu P) 

I -^gdd-^gdd I ][/gl/q2/q3(l + /g4)(l ~ /q5)(l ~ /qd) ~ /gd/qs/qd 
X (l + /gi)(l-/q2)(l-/q3)]}, 

(1) 

and the transport equation for up-quark matter is 

^ _1 f (fp2 (Ppz (Ppi .j 

dt ^ 2EiJ {2n)^2E2 i2TT)^2E3 {2Tr)^2Ei^ ’ 

xS'^ipi +P 2 -P 3 -P 4 ) {go I Mug^ug p [/ql/g 2 (l - /q3)(l + /gd) 

—/q3/g4(l — /ql)(l + /g 2 )] + 5q( 2 I -^UU-S-UU P + | -Adud-s-ud P 
+ I -Aduu^uu P + I ■Ad^(j_j,mj P)[/ql/q 2 (l — /q3)(l ” /qd) 

~/q3/q4(l ~ /ql)(l ~ /q 2 )]} 

1 I' (Pp2 d?P3 d^Pi d^p3 d^pQ 

~2Ei J {2Trf2E2 {2T:f2E3 (2Tif2Ei {2TTf2E5 (2Tif2EQ 

x{2'Kf5^{piEp2 + P3 -Pd -P5 -Pd) I Mugg^ugg P 

X [/ql/g2/g3(l ~ /q4)(l + /g5)(l + /gd) ~ /q4/g5/gd(l ~ /ql)(l + /g 2 ) 
x(l + /g3)] +5q5g(2 I -A^ uug—)-uug P -|- I A^udg-J-udg p + I A^uug-s-uug p 
+ I -A^udg-S-udg I )[/ql/q2/g3(l ~ /q4)(l ~ /q5)(l + /gd) ~ /q4/q5/gd(l ~ fql) 

X (1 — /q2)(l + /gd)] + 9q[-^ I -Aduuu-J-uuu P “*" 4(1 -A^uud-s-uud P 

~t“ I -Adudu—>-udu I ) “t” ^ I -Adudd—J-udd | d”" I -Ad^uu—J-uuu | E~ \ -Ad^^d—>uud I 
+ I -Adudu-s-udu I + I ■Ad^dd->■udd I ”^^1 -Aduuu-^uuu | “^2^1 -A^uud-s-uud I 

+ I -A^udu-J-udu P) + ^ I -A^udd-S-udd H [/ql/q2/q3(l — /q4)(l — /q5)(l ” /qd) 
-/q4/q5/qd(l ” /ql)(l ” /q2)(l ” /qd)]} , 

( 2 ) 

where massless partons have the velocity vector vi and the position vector r; qq and 
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gQ are the spin-color degeneracy factors for the gluon and the quark, respectively; 
the distribution functions for the gluon, the up quark, the down quark, the up anti¬ 
quark, and the down antiquark are denoted by /gi, /ui, /di, fun and /dj, respectively; 
/ui = /di = /ui = /di = /qi- -Adijk-^ijk of which the subscript, i, j, or k, stands for 
a gluon, a quark, or an antiquark is the amplitude for the elastic parton-parton- 
parton scattering. Squared amplitudes were calculated in perturbative QCD with 
different number of Feynman diagrams for different types of scattering—. When 
the squared four-momenta of gluon and/or quark propagators approach zero, the 
divergences of the squared amplitudes for the elastic two-parton scattering and 
the elastic three-parton scattering are encountered. The divergences are removed 
while the propagators are regularized by a screening mass which is evaluated from 
the distribution functions by a formula in Refs.— . It can be proved from the 
transport equations that quark-gluon matter in thermal nonequilibrium eventually 
arrives at thermal equilibrium—. The larger the squared amplitudes, the shorter 
the thermalization time of quark-gluon matter—. 

Every time the elastic two-parton scattering occurs, the two partons’ momenta 
change but the scattering is limited to a plane. Due to the energy-momentum con¬ 
servation only two of the six momentum components of the two final partons are 
free, i.e., angular distributions of the outgoing partons can be defined. The proba¬ 
bility of producing two partons in an incoming parton direction is larger than the 
one in the direction perpendicular to the incoming parton direction. In initially pro¬ 
duced quark-gluon matter more partons are in the gold-beam direction than in the 
direction perpendicular to the beam direction. On the one hand the collision of two 
partons in the beam direction reduces partons in the beam direction and increases 
partons in the direction perpendicular to the beam direction, on the other hand the 
collision of two partons in the direction perpendicular to the beam direction reduces 
partons in the direction perpendicular to the beam direction and increases partons 
in the beam direction. Since more partons are in the beam direction, the scatter¬ 
ing reduces parton number in the beam direction and increases parton number in 
the direction perpendicular to the beam direction. Parton momentum distributions 
gradually become isotropic in such a manner via the elastic two-parton scattering. 

In the elastic three-parton scattering the three initial partons are generally not 
in a plane, so are the three final partons. Every time the elastic three-parton scat¬ 
tering occurs, the three partons’ momenta change. Even restricted by the energy- 
momentum conservation, five of the nine momentum components of the three final 
partons are free. That the parton momentum distributions become isotropic is thus 
expected to be easy via the elastic three-parton scattering. Feynman diagrams ex¬ 
hibit the following processes: a virtual parton breaks into three final partons, two 
final partons, or one final parton and a virtual parton; a virtual parton combines 
with an initial parton to form a final parton; a virtual parton scatters with an initial 
parton to produce three final partons, two final partons, or one final parton and a 
virtual parton; a virtual parton coming from the scattering of two initial partons 
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does not stay in the plane of the two initial partons. Due to these processes the 
three final partons are generally not in a plane. 

Eventually, isotropic parton momentum distributions are established by the 
three-parton scattering and the two-parton scattering. With anisotropic momentum 
distributions of partons in quark-gluon matter created in central Au-|-Au collisions 
at ^snn = 200 GeV at t = 0.2 fm/o^, the transport equations yield isotropic 
momentum distributions of gluons first at t = 0.52 fm/c and of quarks first at 
t = 0.86 fm/c. The thermalization times of gluon matter and quark matter are 0.32 
fm/c and 0.66 fm/c, respectively. The temperature acquired is not close to zero, 
and is much higher than the QCD critical temperature 0.175 GeV. Therefore, the 
early thermalization is a consequence of the two-parton scattering and the three- 
parton scattering. Since the four-parton scattering, the five-parton scattering, etc. 
have certain occurrence probabilities^, they also contribute to the early thermal¬ 
ization. The temperature of the quark-gluon plasma comes from the elastic and 
inelastic two-parton scattering, the three-parton scattering, and the scattering of 
more partons. 
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